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LIST OF SYMBOLS 


axial distance (X /L ) 
height of nozzle (y /L ) 
lateral extent of nozzle (z /L ) 

Jf A 

streamline distance (s /L ) 

distance normal to streamline (n /L ) 

* 2 

mixture enthalpy {h /u m ) 

reference length (throat he'ight - ft) 

Mach number 

« 

th 

molecular weight of i ' species 

h 

mass fraction of i species (frozen mixture) 
number of species (frozen mixture) 
pressure (P'/p^uJ 

JL 

velocity (q'VuJ 

JL ^ 

temperature (T/Tj 

JL 

density (o'Vp^,) 

JL 2 

stagnation enthalpy (H /u m ) 
incremental entropy change 

A 

average po \ ecul ar weight of mixture (frozen mixture - w /wj 

fuel to air equivalence ratio (equilibrium mixture) 

flow inclination 

equilibrium isentropic exponent 

ratio of specific heats (frozen mixture) 

heat capacity at constant volume 

Mach angle 
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SECTION I 
INTRODUCTION 


The process of designing a nozzle for a hypersonic airbreathing vehicle 
involves a complex study of the inter-relationship among many parameters, 
i.e., internal-external expansion, vehicle lift, drag, pitching moments, 
structural and weight limitations. Reference (1) summarizes the various 
stages of the design process, wherein a range of design parameters are 
chosen via a simplified analysis described in Reference (2) with more 
sophisticated calculations performed after the range of parameters has 
been reduced. The second stage of the design procedure (i.e., a more 
sophisticated analysis) is described in Reference (3) and in this report. 

As discussed below the source flow characteristic approach described in 
Reference (3) has been extended and improved. An improved streamline 
interpolation procedure has been incorporated. In addition, all character- 
istic and boundary calculations have been made compatible with frozen, equi- 
librium and ideal gas thermodynamic options, while slip surface calculations 
(cowl interaction) have been extended to underexpanded flow conditions. 

Since viscous forces can significantly influence vehicle forces, pitching 
moments and structura 1 /we ight considerations a local integration via flat 
plate boundary layer skin friction and heat transfer coefficients has been 
included. These effects are calculated using the Spalding and Chi method 
described in Reference (7). All force and moment calculations are per- 
formed via integration of. the local forces acting on the specified vehicle 
wetted areas. 

While this report discusses in detail the improvements mentioned above, 
summaries of appropriate parts of Reference (3) have been included for 
purposes of clarity and continuity. Additional details may be found in 
Reference (3). J 
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SECTION II 
BASIC EQUATIONS 


A. Equations of Motion ~ The equations governing the two 
dimensional, ax (symmetric, or axially expanding inviscid flow of a gas 
mixture in chemical or frozen equilibrium may be written as follows: 

Continuity: (pq) + pq pq pq = 0 

S-Momentum: pq = 0 

dS 3S 

N-Momentunr pq^ ~ + ^ = 0 

3s 3n 

* 

Along a streamline the following relations are valid: 


0 ) 

( 2 ) 

(3) 


Conservation of ^ 12 

Stagnation Enthalpy: ^ = 0 where H = h + j q (A) 


Constancy of Equivalence 
Ratio: 



(equilibrium Mixture) 


(5a) 


or 

Species Conservation: 
and the 

Equation of State: 


3a. 

= 0 (i=1, NSP) (Frozen Mixture) 


p 

p/p = constant • exp (-AS/Cv) 


(5b) 


( 6 ) 


However for a frozen or equilibrium gas mixture (AS/Cv) along streamlines 
is zero in regions of continuous flow. 

J 

B. Thermodynamic State 

Equilibrium Mixture - The static properties are defined via curve fits 

described in Reference (3) i.e., 
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h = h (P ,<t ,T) 
r = r CP,h,<s) 
p = p (P,h,<f) 

Frozen Mixture - The static properties 

Nsp 

h * Z a, h. (T) 
i=1 


Cp f 

If ' C Pf -R/C P _ 


static enthalpy 

(7a) 

1 

equilibrium isentropic 
exponent 

(7b) 

thermal Equation of State 
(density) 

(7c) 

redefined by 


static enthalpy 

4 

(8a) 

ratio of specific heats and 
isentropic exponent (y) 

(8b) 


Nsp 

Cp^: “ S Cp . (T) 


specific heat 


(8c) 


w 


Nsp 

Z 

i=1 


ii 

m. 

i 


-1 


molecular weight 


(8d) 


R 



gas constant 


(8e) 


where pressure, temperature and density are related by 

RJ w 

p =p . — — ~~ thermal Equation of State (9) 

W y M ^ , 

1 co a> J 

the thermodynamic properties Cpj(T), hj(T) are tabulated polynomials, a 
description of which may be found in Reference (4). 
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C, Characteristic System - The system of Equations (1) through 
(6) combined with the thermodynamic system (.7 or 8) yield characteristic 
directions for Mach numbers greater than one. 

Let C denote an up-running and C denote a down-running characteristic. 

+ - 

Then, along a C. characteristic, whose slope is expressed by: 

T 


“ t an (o+p) (lO) 

the compatibility relation may be written 


s I Hu cosy 


d In p t d6+ 


(j £lnl + j 

1 y 2 x 


strip 
cos (0±y) 


dx = 0 (11) 


and along a streamline 


(<&) = tan0 (12) 

» dx 

s 

The local Mach number and Mach angle are given by 

M = q/a (13) 

and ]i = s in 1 ~ (l^) 

where a^ = Tp/p (15) 


It should be noted that a point (x,y) in the flow, properties are completely 
specified by P, H, q,.0, otj or <P. The remaining thermodynamic variables 
are then determined from Equations (7 or 8) 
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SECTION III 

NUMERICAL PROCEDURE AND CHARACTERISTIC NETWORK 


Figure (l) depicts the global grid ordering scheme used in the present 
program. While a free running characteristic network is used the program 
orders and stores data along down-running (C_) characteristics (Reference 
3). The only exception being the initial data line which must be a non- 
characteristic line. The marching proceeds from one C_ line to another 
until the desired flow field is calculated. 

A typical characteristic mesh is depicted in Figure (2), properties being 
known to second order along lines EB, EA and are to be determined at the 
point C. 



The location of streamline point D along characteristic EA or EB represents 
a significant improvement over the procedure described in Reference (3)". 

*A linear interpolation along a characterlst icVs _ second 'order i f the proper- 
ties along the characteristic are computed to second order. Appendix I con- 
tains a proof of the above statement. 


- 6 - 


TR 2! 3 


For Irrotattonal and homentropic flows the streamline location is incidental 
to the calculation of properties at pgint C. However, for highly rotational 
and non-homoen tropic flows the streamline represents a third characteristic 
direction. Thus accurate representation of flow properties along streamline 
DC are necessary for a stable and accurate solution (Reference 5). 



G 1 - 

a tan(0+p) A + 0 tan(0+5) c 

(16a) 


g 2 - 

a:tan(8-p) B + 0 tan(0-p) c 

(16b) 

and 

G 3 = 

a tan 0^+0 tan 0^, 

(16c) 


The Mach angles (p) are the local values corresponding to either frozen or 
equilibrium flow. The a and 0 in the above and following equations are 
used as artifices in averaging properties along characteristics. In a first 
approximation a would be set equal to one and 0 equal to zero. Once proper 
ties at point C are determined, the coefficients involved in the calculation 
are averaged by setting both a and 0 equal to one~half. This corresponds to 
the second iteration. Writing Equation (10) in finite difference form 


V y A 

X C _X A 


and 


V y B 

V X B 


= G, 


Solving the above yields 


X C = 


y B~ y A +G 1 X A" G 2 X B 


G r G 2 


(17a) 


(17b) 


J 

(18a) 


y C = y A + G 1 (X C“ X A } 


(18b) 
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A? = BC s.lnu co«w ) + p C sln H co s y . ) 


( 19 a) 


0 - ^sinu cosu^ + p ^slnu cosyi^ 

B Y i 


( 19 b) 


The Mach number and angles u are the local equilibrium or frozen values. For 
equilibrium flow the ratio of specific heats y is replaced by the equilibrium 
exponent r. Similar remarks apply to all further characteristic coefficients. 


„ , r /Sin0 sinu \ 

A 2 ■ J 1 + 

A 

„/sin0 sinu \ t 
y cos(0+u) 

C 


(20a) 

r /COS0 sinu \ 

+ J 2 ta( xToV(W) ) + 

A 

„/cos0 sinu \ i 
P^x cos (0+u) J 
C 


(20b) 

B 2 J 1 “ y cos (0-y) + 

B 

/ 5 in 0 sinu \ ■, 
V cos (0 - u) J 

C 


(20c) 

. , r /cos0 sinu \ 

+ J 2 l “ ( x cos (e-p) + 

,cos0 sinu \ i 
'x cos (0-y) J 


(20d) 

B 

C 



Hence, along AC (C + characteristic) 




A 1 (ln P C ‘ 1n P A } + 0 C" 0 A 

* A 2 W * 

o- 

(21a) 

and along BC (_C_ characteristic) 


1 


B 1 (ln P C " ln P B } ' V 6 B 

+ B 2 (x c -x b ) = 

s 

0 

(21b) 
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Solving the above equations for In yields 

In P C = (A, ,n P A + B 1 ln P B V 0 B “ (A 2 +B 2 ) x C 

( 22 ) 

+ A 2 X A + 8 2 X B^ A 1 + B 1^ 
hence p c = exp (In p c ) and the flow Inclination is 

0 C = e A - A, (In p c - In p A ) - A 2 (x c -x fl ) ( 23 ) 


if either the entropy or stagnation enthalpy Is not constant throughout 
the flow field, the streamlines are the third family of characteristics, 
whose slope is given by 


4* = tan 0 


in difference form, 


( 12 ) 


y C~ y D 

V X D 


= G. 


(24a) 


Referring to Figure (2), the point D lies between points EA or EB and can be 
located by an iterative procedure using Equation (24a). Properties at point 
D are then obtained by linearly interpolating along the appropriate charac- 

p 

teristics EA or DB, location of point D yields the variables Pj , 0J , (P/p ) d 
Pd» 9d» $ d or ctj d (M ,NSP) . 


From the streamlint Ccis^-vation Equations (4), (5a) or (5b) and the Equation 
of State (6) 

H C = h d 

+ C = 

or a tr * a , (1 = 1, NSP) 
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and 


P r =(P c (/P/p r ) ) 
C D 



Note that as a first approximation T c |v Tp. 

If the flow is in frozen equilibrium the remaining properties at C are 
determined directly from 


W = 
c 


r NsP 

2 


T "I 


L M 


R - 2° 
c W,- 


T = 
c 


A 

v p ' 


V K 2 

CO 00 

w 


C P - 


NsP 

s «*<:* C Pi (T c> 

H -4 

t =1 


NsP 

h c = 2 a * c * (T,») 

i=1 


[jfc(H - h ) ] 


Y c = 


'Pc 


C - — 

°P C Cr 


If the flow is in chemical equilibrium an iteration is required to determine 
the remaining properties at C since the temperature i s obtained by inverting 
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the ^urve fit (7’a) l.e., h = h (P,4>,T) consistent with the density obtained 
from the equation of state. 

Knowing the temperature »' then yields the static enthalpy h . The velocity 

w C 

q Js obtained from 
c 

q c - [2 (H c - h c )J 1 

and the curve fit (,7b) yields 

r c = r (p,h,<s) 

Equations (13), (14) and (15) then yield H c >PQ> a c - 

The complete calculation ?s then repeated with all coefficients averaged, by 
setting a and 0 equal to 1/2. if properties change significantly between 
these two sets of calculations, this generally implies too large a mesh in 
this region. 
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SECTION IV 

BOUNDARY CALCULATIONS ■ 


The following boundary calculations are permitted in Program NOZBOD; addi- 
tional details may be found in Reference (3 ) • 

A. Upper or Lower Wall - The nozzle wall shapes, either upper wall 
(cowl) or lower wall (vehicle undersurface) are specified by ploynomials of 
the form 

*, 2 

yj = Aj (x-x.) + Bj (x-x.) + Cj ( 25 a) 

where x. is the origin of the wall. A maximum of 5 wall segments are permit- 
ted, i.e., j = 5. The wall slope 0 w i ? given by 

tan 0 = 2A. (x-x.) + B. (25b) 

w J i J 

In Figure ( 3 ), DC is the specified upper nozzle wall. 



FIGURE 3. WALL POINT 

Given the flow deflection ( 0 ) at point C the pressure P c is computed from 
the characteristic relation along AC (C+ char). 
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B. Shock Phenomena - The program developed has the capability of 
computing the shock strength associated with an inviscid supersonic over 
or under-expansion process, and a shock propagating into a nonuriform media. 


1. Hugoniot Relations - Assume a coordinate system oriented 
along (t direction) and normal to the shock surface (n direction) as shown 
in Figure (4) . The angle sigmao is the direction cosine of the shock with 

r\j 

respect to the Cartesian direction x, and u and v are the velocity components 
in the n and t directions. 


n 

- -sinai + 

X 

cos o i 

y 


t 

« COSO i + 

X 

s i no i y 

-► 

\> * % 

A A 

V 

* u r. + v f 

t * ui + 

X 


(26a) 

(26b) 

(27) 



FIGURE 4. SHOCK COORDINATES 


The Rankin Hugoniot relations for a mixture in chemical or frozen equilibrium 
take the form 
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Cont Inui ty : 



P 2 U 2 


( 28 ) 


Normal 

Momentum: 


>v2 

P, ♦ P,u, 


<v2 

P 2 + P 2 U 2 


(29) 


Tangential ^ y 

Vf ■ V 

Momentum: tj 


(30) 


Energy: 


H - h + 


2 


constant 


(3D 


State: 


p = p (p , h , 4> 1 ) 


P 


Pw 

R T 
o 


w 

OO 



(Equilibrium Mixture) 
(Frozen Mixture) 


(32a) 

(32b) 


where 


<t> = constant 


(Equilibrium Mixture) 


(33a) 


u. * constant 


(Frozen Mixture) 


(33b) 


Employing the jump relations for a given shock angle and upstream conditions 
requires an iteration process since the mixture is calorically imperfect. 


Let 1 designate upstresm conditions and 2 downstream conditions. To solve 

«v 

the jump relations knowing conditions at 1 , a value for u 2 is assumed. The 
density p 2 is computed using Equation (28), p 2 is computed using Equation 
(29) and Equation (31) yields a value for h 2 - The State Equation (32) then 
yields an alternate value for the density. If this value for density does 
not agree with that calculated from continuity to within a specified toler- 
ance, the value of u 2 is perturbed and this process is repeated until con- 
vergence is achieved. 
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2. Shock Point Calculation - Referring to Figure (5 ) » a 
typical shock wave calculation is performed as follows. A value of the shock 
angle o ^ is assumed, and a simultaneous solution of the equations 

. B 



and 


yields 


V y A 


= \ (tan c c + tan a ) 


C A 


V y B 

X C _X B 


Since this flow is nonuniform a characteristic calculation simTlar to an in- 
terior point calculation yields the flow properties at Cj. Note that point 
E on the C + characteristic EC^ is interpolated along FA^ . The jump relations 
(Equations 28 - 33) are solved using the determined upstream conditions based 
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on the assured angle o^.. This yields all properties at Using the 

deflection angle 0- calculated from the jump conditions, a C character- 

L 2 + 
istic calculation performed along (G-Cj) yields an alternate value of the 

pressure at p^- The pressures are compared anu if the difference exceeds 
a specified tolerance, the value of c^. is perturbed and the process repeat- 
ed until convergence is obtained. After convergence with a ■ 1, 6=0, 
the complete calculation is repeated using a = 1/2, B = 1/2. 


3. Shock Reflection at Wall - The incident and reflected strength 
of a shock wave at a wall boundary is determined by the condition that 
downstream of the reflected wave (3) the flow deflection at the wall must 
equal the wall slope; Fic,ure (6) depicts this interaction. 



C. Cowl Interactions 


1. Under-expansion Interact ' on - The program developed has 
the capability of computing the under-expansion interaction produced by pres- 
sure mismatch between the nozzle and a surrounding airstream. This situation 
is depicted in Figure (7a). Under-expansion conditions occur as a result of 

P. > P or 0. > 0 or some combination. Generally P. >> P defines an under- 
j e j e J e 

expanded flow. 
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2. Over-expansion Phenomena - The nozzle over-expansion at 
the cowl is computed in a similar fashion to the under-expanded phenomena, 
except that a shock wave is required in the nozzle flow and an expansion in 
the external flow as depicted in Figure (7b). 




FIGURE 7b. OVER-EXPANSION INTERACTION 
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The jet flows are assumed to be H^'air combustion products either in chemi- 
cal 01 frozen equilibrium. The external flow is assumed to be a perfect 
gas . 


0. Contact Surface - A contact surface is a stream surface of 
the flow, therefore, the pressure and flow deflection must be continuous 
across the discontinuity. Figure (8) illustrates a contact surface calcula 
tion for supersonic flow. In the present program the characteristic on the 



FIGURE 8. CONTACT CALCULATION 


external flow side is replaced by a Prandt 1 -Meyer pressure-flow deflection 
relation or shock expansion theory and the external flow is assumed uniform. 
The solution requires a iterative procedure similar to a wall boundary calcu- 
lation except that the shape of the boundary (i.e., contact surface) is not 
* 

known a priori . 


i 
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* SECTION V 

THRUST, LIFT AND PITCHING MOMENT 

The following definitions are used in this report for thrust, lift, pitching 
moment 


T 


L 




vis 


I 

y 


dA 


"vis 


(3M 


(35) 


M y 



(36) 


Figure (9) gives the orientation of the vectors with respect to the vehicle. 
Internally the integrals range over all the vehicle surface areas. Ex- 
ternally they range over the complete vehicle undersurface as defined by 
the bounding stream surface and/or flow fence. 



FIGURE 9- THRUST, LIFT, MOMENT 
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As defined above the nozzle 
forces. These effects will 

as 


thrust, lift and pitching may include viscous 

be discussed below. However, we define T , , L , 

’ vis’ vis 


T 


vis 



(37) 


L 


vis 



(38) 


Surface Area Computation - Internally the nozzle surfaces are 
defined as the cowl internal surface, vehicle undersurface, and the nozzle 
sidewall. Externally the surfaces are defined as the nozzle undersurface 
and/or flow fences. The lateral extent of the cowl surface and nozzle under- 
surface is defined by the degree of lateral expansion desired. 

It is assumed that the nozzle area can be approximated by a series of ele- 
mental quadrilaterals, as shown in Figure (10). From Reference (6) the unit 
normal for an elemental area may be obtained by defining two surface tangent 
vectors from the diaganols of the quad i 1 at era 1 . 



FIGURE 10. ELEMENTAL SURFACE AREA 
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That is 


T . ■ T 1x * 

T, I 
iy y 

* T 1Z '* 

(39«) 

^ 2 * T 2y * 

* 

T 2y ‘y 

♦ t i 

2z z 

(39b) 

where 




T 1x • *3 - '2 ■ 

T iy 

*3 ' V T lz * z 3 ’ Z 1 


T 2x * \ * V 

T 2y 

* V T 2z ‘ z l. " z 2 


and the normal N is defined 

as 



N - T 

2 xT 1 



and the unit norma) as 




A 

-► 

N 


m 

n ■ 

1 i 1 



A tangent plane is constructed using the normal vector and the two 

tangent 

vectors T^, . The corners of the 

surface element are projected onto chis 

plane anH the area and centroid of 

the quad! lateral are calculated 

as described 

In Reference (6) . 


• 


Typical nozzle elemental areas are 

shown in Figures (11) and (12) 
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FIGURE 11. UPPER OR LOWER SURFACE AREA ELEMENTS 




Lower 

Suface 




FIGURE 12. SIDEWALL AREA ELEMENTS 
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SECTION VI 
VISCOUS EFFECTS 


Local skin friction and heat transfer coefficients are computed via curve 
fit data supplied from Reference (6). These fits are based on the Spalding 
and Chi method of Reference (7). That is a suitably transformed skin-friction 
coefficient is given by incompressible formulas based on a suitably trans- 
formed Reynolds number, i.e.: 


where 


C - C /F 

T « T l C 


C f ^ - f (Rxj) , Rx, - F^ • Rx 


0 * 1 ) 

0 * 2 ) 


■ local skin friction coefficient 


Rx ■ Reynolds number 


( )j » indicates incompressible 
( )^ ■ indicates compressible 


Now for Rx. > 25^*0 the local skin friction is given from Reference (8) as 

3 

C fj - .088 (log Rx j - 2 . 3686) / (log Rx-. - 1.5) ( 1 *3) 


and frou reference (5) 

F c - A/ [ ARSIN (^-) ♦ ARSIN (^) ) 


(V.) 


h aw 

A - - 1 


0*5a) 
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The local properties external to the boundary layer are the local data com- 
puted by NOZBOD and are assumed to act through the centroid of the eleme 
area computed above. The computation requires that a boundary layer oric , 
be specified since the nozzle is assumed to be an extension of the combustor. 
In addition a recovery factor for an adiabatic wall calculation is required. 
However, as a user option wall temperature distributions may be specified. 

Local heat transfer coefficients are computed from a modified Reynolds 
analogy for turbulent flow 


St - Sh • C f /2 

The program requires "Sh" as an input item. 


(*7) 


J 
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SECTION VI I 

CHARACTERISTIC COALESCENCE 


Since the program developed stores data on c characteristics, shock wave de- 
tection for both first and second family shocks i; straight forward. A shock 
wave is present if two characteristics of the same family intersect each other 
within the flow field. 

While characteristic crossings are detected automatically the program does not 
insert shock waves into the flow field calculation. However, by judiciously 
eliminating data points and noting these locations an approximate shock wave 
trace may be obtained during the course of the calculation. 

For down-running shocks, detection must take place prior to the calculation 
of a new data point. Figure (13) illustrates the procedure. Since the flow 



i 

FIGURE 13. DOWN-RUNNING (c‘) SHOCK DETECTION 
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field is known on Line 1, and up to point K on Line J, and if point K+1 is 
the next data point to be calculated, a crossing of c characteristics may be 
determined from the following relations. 


<vV 


"•<vV 


where 


and 


(y c“ y m ) 


N* (x -x ) 
c m 


1 


M 


31 



k-1 

k-1 


y-y m _, 

m m- i 

x -x - 1 
m m 


If x c >x j^ then the initial shock strength is related to the relative distance 
(S^-S^AS^) where S is the running distance along the c characteristic and 
AS^ is an average n.osh distance on Line 1. If the above ratio is less than a 
given tolerance no further data on Line J is calculated and the new initial 
data line becomes Line J up to point K and the remainder of Line 1 from m+1 
to m | a5t ’ A new data line is calculated from this hybrid data and the process 
repeated . 


Detection of up-running shocks may be accomplished a^ter completing the calcula- 
tion of new data Line J. Figure (14) is an illustration of the process. The 
crossing of c + characteristics may be obtained from the following relations 


<V C -Y k > - "-(x c -y k ) 


( VW ’ N ' (x c' x k*l ) 


where 


y k+r y m +1 

X k+r X m-* 1 
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and 


N 


V y 

x. -x 
k 


m 

m 


The up-running (c + ) initial shock strength is determined by the relative 

distance (x -x, )/(y ,-x ). If the above ratio is less than a given tolerance, 
c k ''km 

the point k+1 is dropped from the data Line J and the next data line is cal- 
culated using Line J minus point k+1. 


Thus, a continuous face of both up-running and down-running shock waves may be 
obtained. It should be noted that the process of dropping selective data 
points is valid only for shock waves which do not appreciably alter the flow 
field entropy distribution. 
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SECTION VIII 
SUMMARY AND CONCLUSION 


In Reference (3) a unique characteristic procedure was developed which 
computes realistic scramjet nozzle exhaust flow fields, using an axially 
expanding coordinate system (line source) allowing lateral nozzle area 
variations to be accounted for in a quasi two dimensional fashion. Addi- 
tional geometric flexibility was incorporated in the numerical procedure 
through the use of additional coordinate systems i.e., axisymmetric and 
Cartesian. The technique is second order in a characteristic sense. 

The present study reported herein extends and improves upon the authors 
previous work. Improved streamline interpolation procedures insure nu- 
merical stability and improve program accuracy. The extension of the 
present work to include plume shape affects on vehicle forces removes a 
major shortcoming of the previous work. In addition, sophisticated force 
and moment calculations allows for accurate incorporation of s idewal 1 /fence 
forces and viscous effects (internal drag and heat transfer). Thus heat 
transfer effects can be used to assess vehicle structural and weight con- 
siderations. Shock detection techniques have been incorporated so as to 
accurately determine where body shapes should be modified in order to 
allow shock free flow. 

Thus, it is felt that the current program will give the designer of scram- 
jet nozzle exhaust flow fields increased flexibility not previously avail- 
able. 


) 
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APPENDIX 

Linear Interpolation on Characteristics 
A Second Order Procedure 


.'b 

•J V 

IV 

'*ii 


Employing a 1 i near interpolation procedure ona characteristic line calcula- 
ted to second order jj_ consistent with a second order algorithm. This point 
can be inferred from Ferri's article in Vol . 6 of the Princeton series but 
apparently, is not universally accepted (i.e., Reference 5). Hence, a simple 
proof' of this statement is presented. 

Along any line AC, a series expansion for the pressure and flow deflection 

* 

are wri tten: 


ilii 


P C - P A + (P x> A 4x AC + < P xx\ ^ /2 \ c * 0 <4 * 3) AC 


(ai ) 


hi 


e r « + (0 V ) Ax + (e ) (Ax 2 /2) + 0 (Ax 3 ) 

L M X A AU XX A AC AC 


(A2) 


but 


(P ) = (P x ) + (P vx ) AX. C + 0 (Ax*) 

X C X A XX A A AC 


(A3) 


(e ) = ( 0 ) + (e vx ) Ax.. + 0 (Ax J 

X c A xx A AC AC 


m 


(where x denotes distance along AC). 


,00. 



Substituting (3) into (1) and (b) into (2) we obtain 

P C - P A + I(P X> A + ' P xh c (Ax/2) ac - o (Ax 3 ) ac , 

8 C - e A + t(8 x ) + (e x )] (ax/2) ac + 0 (ax 3 ) ac 

A C 


(A5) 

(A6) 


A-1 


I 
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The above expressions are valid along any line AC. Nov/ at points A and C, 
we write the compatibility relations, assuming that AC is a downrunning char- 
acteristic. 


**\ - A A CP x\ = ° 
" A C (p x } c = ° 


(A7) 

(A8) 


Solving the system (5) - (8) for (P ) + (P ) , we obtain 


*A X C 


<P x' A - 5 T (A c -A/ W - A C (P <f P A> 


(A9) 


and 


(PJ 




{0 c" e A ) “ A A (P C- P A } 


(A10) 


Now consider a point x’ between A and C. The pressure at this point, to 
second order is given by 

^ + (Px) r -(Px). -v 2 9 

' C -A (x * -x) +Q {Ax 3 } 


P" = P. + (P v ) (x* -x) + 

* A A 


Ax 


(All) 


AC 


where (P ) and (P ) ere given by (9) and (10). Till this point, all rela- 
X A X 

tions are quite general and have not required that a second order compatibility 
relation exist between A and C. We now make use of this relation by stating 
that by second order, we imply that the relation 

A +A 

( 0 C - e A ) - (P c - P A ) = o (A12) 

) 

is satisfied between points A and C in a convergent fashion as detailed in 
Reference (9).- Then, substituting Equation (12) into Equations (9) and (10) 
we obtain: 


A-2 
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(A1 3 ) 
(A14) 

Hence, substituting these relations into Equation (11) yields 

* (X*-X ) 

P * P A + (P C “ P A } 7 VX7 + 0 ^ 


(P x J , 


(P X } 


and 




C p c -P A ) 

W 


which clearly demonstrates that a linear interpolation for pressure (or 
flow deflection) on a characteristic calculated to second order is con- 
sistent with a fully second order approach. 


j 
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